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Observational evidence for a connection between supermassive 
black holes and dark matter haloes 
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ABSTRACT 

We present new velocity dispersion measurements of sample of 12 spiral galaxies for which 
extended rotation curves are available. These data are used to refine a recently discovered 
correlation between the circular velocity and the central velocity dispersion of spiral galaxies. 
We find a slightly steeper slope for our larger sample, we confirm the negligible intrinsic 
scatter on this correlation, and we find a striking agreement with a corresponding relation 
for elliptical galaxies. We combine this correlation with the well-known M BH -<r relation 
to obtain a tight correlation between the circular velocities of galaxies and the masses of 
the supermassive black holes they host. This correlation is the observational evidence for an 
intimate link between dark matter haloes and supermassive black holes. Apart from being an 
important ingredient for theoretical models of galaxy formation and evolution, the relation 
between Mbh and circular velocity can serve as a practical tool to estimate black hole masses 
in spiral galaxies. 
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1 INTRODUCTION 

Soon after the discovery of quasars in the 1960s, it was recog- 
nized that only supermassive black holes (SMBHs) in the centres of 
galaxies could provide the mechanism to feed active nuclei. An in- 
triguing problem arose however, when it became clear that quasars 
were much more abundant in the early universe: what happened to 
these SMBHs ? A natural explanation is that SMBHs still reside 
in the centres of present day galaxies, but somehow the fueling of 
the active nucleus stopped, for example due to the competition with 
star formation (Di Mateo et al. 2003). The detection of SMBHs in 
quiescent galaxies is a difficult task. Except for some special cases 
such as NGC 4258 (Miyoshi et al. 1995) and the Milky Way (Ghez 
et al. 1998; Genzel et al. 2000), the detection of a SMBH depends 
on stellar or gas kinematics. In order to be able to actually detect 
the SMBH, it is necessary to resolve the sphere of influence, i.e. the 
central region where the black hole starts to dominate the potential 
of the galaxy. With the limited resolution of ground-based spectro- 
graphs, this was possible for only a handful of nearby galaxies. This 
changed drastically with the advent of the HST, whose excellent 
spatial resolution enabled to resolve the SMBH for a few dozen of 
nearby galaxies (see Tremaine et al. 2002 and references therein). 
In nearly all of the investigated galaxies, SMBHs were discovered 
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with masses roughly between 10 7 and 10 9 M©, Very recently, the 
detections of SMBHs of much lower masses have been reported in 
globular clusters (Gebhardt, Rich & Ho 2002; Gerssen et al. 2002, 
2003), although these detections are still controversial (Baumgardt 
et al. 2003a,b; Ho, Terashima & Ulvestad 2003). 

Having detected and measured SMBHs for a sizable sample of 
nearby galaxies, we can proceed to tackle more fundamental ques- 
tions concerning their formation and evolution. A convenient way 
to do so is by studying the relation between SMBHs and the galax- 
ies that host them. It was found that black hole masses are corre- 
lated with parameters of the bulges of their host galaxies. Nearly 
ten year ago, Kormendy & Richstone (1995) found a correlation 
between the mass A/bh of the SMBH and the blue magnitude Lb of 
the bulge. More recently, a tighter correlation between Mbh and the 
velocity dispersion a of the bulge was discovered independently by 
Gebhardt et al. (2000) and Ferrarese & Merritt (2000). Finally, Gra- 
ham et al. (2001) presented an intriguing correlation between Mbh 
and the so-called concentration parameter of the bulge (or Sersic 
index), in the sense that bulges with more massive black holes have 
steeper cusps. 

This apparently tight link between bulges and SMBHs reflects 
an important ingredient that should be reproduced (and hopefully 
explained) by theoretical models of galaxy formation. Actually, the 
tightness of the correlations mentioned above is somewhat sur- 
prising. Indeed, in most of the state-of-the-art models (e.g. Silk 
& Rees 1998; Adams et al. 2000; Kauffman & Haehnelt 2000; 
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Burkert & Silk 2001; MacMillan & Hendriksen 2002; Wyithe 
& Loeb 2002; Volonteri, Haardt & Madau 2003; Di Matteo et 
al. 2003), the total galaxy mass (or dark matter mass A/dm), rather 
than the bulge mass, plays a fundamental role in shaping the 
SMBHs. A close correlation could therefore be expected between 
Mbh and Mdm, rather than between A/bh and the properties of the 
bulge. 

Unfortunately, this hypothesis is not straightforward to test 
observationally, because the measurement of A/dm and A/bh are 
difficult. An interesting approach to tackle this problem was taken 
by Ferrarese (2002). She argued that a correlation between A/bh 
and A/dm should be reflected in a A/bh - v c correlation, where v c is 
the circular velocity in the flat part of the rotation curve of spiral 
galaxies. Unfortunately, there are only four spiral galaxies with se- 
cure SMBH masses (Genzel et al. 2000; Miyoshi et al. 1995; Bacon 
et al. 2001; Lodato & Bertin 2003), and only two of them have a 
measured rotation curve. Therefore, Ferrarese used the tight A/bh - 
a correlation to estimate black hole masses in a larger sample of 
galaxies: in fact, from a literature search, she found a tight corre- 
lation between v c and a (with both v c and a measured in km s~ ), 

logv, = (0.84 ± 0.09) log<r + (0.55 ±0.19). (1) 

The tightness of this correlation (it has a negligible intrinsic scatter) 
suggests an intimate link between SMBHs and dark matter haloes. 
A possible caveat in Ferrarese's analysis is that it is based on only 
16 spiral galaxies, including three spirals with a < 70 kms - 
which do not satisfy the relation Q. A larger database of galax- 
ies with extended rotation curves and velocity dispersion measure- 
ments would be useful to check and refine this correlation. This 
is the goal of the present Letter. We discuss the sample selection, 
the observations and data reduction in Section 2. We analyze and 
discuss these data in Section 3. 



2 THE DATA SET 

We aim to increase the number of galaxies with accurate measure- 
ments of the circular velocity and the central velocity dispersion, 
in order to investigate a possible link between SMBHs and dark 
haloes. It is a well-known fact that the velocity dispersion of galax- 
ies is a fairly difficult quantity to measure accurately: values for the 
velocity dispersion quoted in the literature often differ consider- 
ably depending on the authors' favorite fitting techniques, template 
stars, absorption line, signal-to-noise, etc. In order to construct a re- 
liable sample of galaxies, we therefore chose to select galaxies from 
a homogeneous sample with measured extended rotation curves, 
and to measure the velocity dispersions in a consistent way. 

Our galaxies were drawn from the sample of Palunas & 
Williams (2000), who present a sample of 74 southern spiral galax- 
ies. For each of these galaxies they constructed axisymmetric mass 
models based on accurate two-dimensional Ha velocity fields and 
/-band imaging. 32 of these galaxies have rotation curves ex- 
tending to beyond the optical radius , 26 of which are reason- 
able smooth. These formed an ideal target for our study. We were 
granted one night of observing time (4 May 2002) to measure the 
velocity dispersion of these 26 galaxies at the ESO 3.6m telescope 
in La Silla. Unfortunately, due to a very poor sky transmission, the 

The optical radius used by Palunas & Williams (2000) is -R23.5, which 
is the radius of the isophote corresponding to an /-band surface brightness 
of 23.5 mag arcsec -2 
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NGC 2844 
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NGC 2903 


SBbc 


180 ±4 


106 ±13 


7.02 ±0.24 


NGC 2998 


SBc 


198 ±5 


113±30 


7.13 ±0.48 


NGC 3198 


SBc 


150 ±3 


69 ±13 


6.27 ±0.36 


NGC 4062 


SBc 


154 ±13 


90 ±7 


6.73 ±0.19 


NGC 4258 


SBbc 


210 ±20 


138 ± 18 


7.48 ± 0.24 


NGC 4565 


Sb 


264 ±8 


151 ± 13 


7.64 ±0.17 


NGC 5033 


Sc 


195 ±5 


122 ±9 


7.27 ±0.16 


NGC 6503 


Sc 


116±2 


48 ±10 


5.64 ±0.42 


NGC 7331 


Sbc 


239 ±5 


139 ± 14 


7.50 ±0.19 



Table 1. The galaxies in our sample. The galaxies above the horizontal 
line are the galaxies from Palunas & Williams (2000) for which we have 
obtained a velocity dispersion. The galaxies below the horizontal line are 
the galaxies from Ferrarese (2002) with a rotation curve extending beyond 
the optical radius. The first two columns contain the name and morpholog- 
ical type of the galaxies. The third and fourth columns contain the circular 
velocity and bulge velocity dispersion of the galaxies, with errors. The last 
column contains an estimate of the SMBH mass, based on the Mbh-"" 
relation (see text). 

exposure times had to be multiplied by a factor 2.5 to obtain a suf- 
ficient signal-to-noise, such that we could only observe 12 galaxies 
from the sample. The observed galaxies are listed in Table 1. 

We have taken long-slit spectra of these 12 galaxies with the 
EFOSC2 instrument through a 0.5 arcsec slit, giving a spectral res- 
olution of 0.31 nm FWHM. We used the Gr#08 grating, which en- 
abled to measure the velocity dispersions of the galaxies from the 
Mgl and Fe lines around 520 nm. The exposure time varied from 
25 min to 150 min for the faintest galaxies. Standard data reduc- 
tion and calibration of the spectra were done with the ESO-MIDAS 
package. This includes bias subtraction and trimming, cosmics re- 
moval, wavelength calibration (Legendre order 4 is used), removal 
of the S-distortion/tilt of the slit, sky subtraction and airmass- 
correction. No effort was done to flux-calibrate the spectra, as we 
are only interested in the velocity dispersion. The reduced spectra 
have a wavelength domain between 430 and 630 nm, with a pixel 
scale of 9.8 nm/pixel. 

The velocity dispersion a of the galaxies was determined by 
means of a direct % 2 fitting technique, which fits line profiles di- 
rectly to the spectra after convolution with a stellar velocity tem- 
plate spectrum. The fitting routine calculates the full covariance 
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Figure 1. The correlation between the circular velocity v c and the central velocity dispersion <r. The left plot shows the v c — a correlation for the 24 spiral 
galaxies with rotation curve beyond the optical radius and a velocity dispersion <r > 80 kms -1 . The green dots are the data points from the sample of 
Ferrarese (2002), the red dots are the galaxies with new velocity dispersion from this paper. The right plot zooms out and adds the elliptical galaxies from 
Kronawitter et al. (2000) to this plot. The straight line in both panels represents equation 0. 



matrix, from which we could determine the error bars on the disper- 
sion. The spatial aperture over which the fit was done was chosen 
to maximize the signal-to-noise ratio of the spectrum, which was 
targeted to be of the order of 15. In fact, our signal-to-noise values 
ranged between 11 and 35. We found that the value of the disper- 
sion did not depend sensitively on the extent of the aperture. The 
velocity dispersions of our 12 galaxies can be found in the upper 
part of Table 1 . 

The circular velocities of galaxies were taken from Table 1 
in Palunas & Williams (2000). These are determined by taking a 
weighted average of the rotation curve points where the rotation 
curve becomes flat. We repeated this exercise on the original data, 
which were kindly provided in tabular form by P. Palunas, to get 
an estimate for the errors on the circular velocity. We find identical 
values as Palunas & Williams (2000), with errors of the order of 5 
to 10 per cent. 



3 ANALYSIS AND DISCUSSION 
3.1 The correlation between v c and a 

Combining our 12 galaxies with the 16 spiral galaxies from the 
compilation of Ferrarese (2002) with a rotation curve extending 
beyond the optical radius (see Table 1), we obtain a total sam- 
ple of 28 galaxies. In the left panel of figure Q we plot the cir- 
cular velocity versus the velocity dispersion for these 28 galaxies. 
For the 24 galaxies with a velocity dispersion greater than about 
80 kms -1 , there is a very tight correlation between v c and a. We 
fitted a straight line to these data, taking into account the errors on 
both quantities (Press et al. 2002), and obtained 

\o%v„ = (0.96 ±0.11) log u+ (0.32 ±0.25). (2) 

In this correlation, both v c and a are expressed in km s ~ 1 . However, 
in a linear regression analysis, this is not the most meaningful unit 
to express these quantities. Tremaine et al. (2002) argued that the 
choice of appropriate units is important to avoid a strong correlation 
of the errors in slope and zero-point. When we repeat our fitting 



routine after choosing a new unit uq = 200 km s 1 , we obtain as a 
best fit 

log f — ) = (0.96 ±0.11) log ( — ) +(0.21 ± 0.023), (3) 
\ / \ Mo J 

where indeed the uncertainty on the zero-point is much smaller. 

Our newly derived v c - a relation has a slightly larger slope 
than the original correlation Q found by Ferrarese (2002) based 
on 13 galaxies. Both slope and zero-point are consistent within the 
one sigma error bar, though. The tightness of the correlation is still 
astonishing: we find a reduced \r = 0.281, corresponding to a 
goodness-of-fit of 99.9 per cent. The v - a relation can hence be 
regarded as having a negligible intrinsic scatter. Important, more- 
over, is the fact that this correlation appears to be robust: even with 
the sample of galaxies nearly doubled, there are still no significant 
outliers. Ferrarese's v c - a correlation broke down for galaxies with 
dispersions below about 80 kms -1 . We confirm this result: the only 
galaxy we observed with a < 80 km s~ (ESO 383-G88) joins the 
three galaxies from Ferrarese's original sample (M33, NGC3198, 
NGC 6503 and ESO 383-G88) in lying significantly above this cor- 
relation. 

It is interesting to compare this correlation to a correspond- 
ing one recently found for elliptical galaxies. Based on stellar dy- 
namical models for 20 ellipticals constructed by Kronawitter et 
al. (2000), Gerhard et al. (2001) discovered a very tight relation 
between the central dispersion and the circular velocity (the circu- 
lar velocity curves of ellipticals were found to be flat to within 10 
per cent). Ferrarese (2002) reports for these 20 ellipticals the cor- 
relation 

log« c = (0.94 ±0.11) log a + (0.31 ±0.26), (4) 

with a reduced Xr = 0.66. Both the slope and zero-point of this 
correlation agree amazingly well with the v c - a correlation of 
our spiral galaxy sample. These expressions represents a nearly di- 
rect proportionality between the bulge velocity dispersion and the 
halo circular velocity, with the proportionality constant about two 
thirds. To first order, elliptical galaxies form a dynamically uni- 
form family, such that we expect a proportionality between a and 
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Figure 2. The correlation between the circular velocity v c and the SMBH mass Mbh- Black hole masses are estimated from the velocity dispersions through 
the Mbh - a relation, as given by equation J5J. The layout is similar to figure^ The straight line in both panels represents equation 



v c (Gerhard et al. 2001). For spiral galaxies, however, this propor- 
tionality cannot be explained by simple dynamical arguments, as 
convincingly argued by Ferrarese (2002). The fact that there is such 
a strong correlation between both velocity scales in spiral galaxies 
is hence not obvious, and indicates a fundamental correlation in the 
structure of spirals. Moreover, the fact that both spiral and elliptical 
galaxies obey exactly the same correlation is absolutely striking. 

3.2 The correlation between A/bh and v c 

Although derived from mainly late-type spirals (all galaxies in our 
sample except NGC 2844 have a Hubble type of Sb or later) with a 
between 90 and 180 km s~ , the v c -a correlation appears to hold 
as well for a larger class of galaxies, ranging from ellipticals to late- 
type spirals, with a dispersion range up to 350 kms -1 . A similar 
situation appears to apply for the A/bh —o relation. From the three 
recently discovered empirical relations linking A/bh to the bulge 
properties of the host galaxies (see Introduction), the A/bh - a rela- 
tion is the tightest one. The two original papers describing this cor- 
relation found significantly different slopes: 4.8 ± 0.5 (Ferrarese & 
Merritt 2000) versus 3.75 ±0.3 (Gebhardt et al. 2000). There has 
been a vivid discussion since over the exact slope of this relation, 
and it is argued that these differences can be ascribed to different 
fitting techniques, different samples and different measures of the 
central velocity dispersion (Merritt & Ferrarese 2001a, b; Tremaine 
et al. 2002). The most up-to-date version of the A/bh-c correla- 
tion, from Tremaine et al. (2002), reads 

log = (4-02 ± 0.32) log (j-j + (8.13 ± 0.06) (5) 

This relation is derived from empirical data for some 30 nearby 
galaxies galaxies with secure SMBH detections. Most of these 
galaxies are elliptical galaxies with only a minor fraction of spi- 
ral galaxies present. Nevertheless, this correlation seems to hold 
for lenticular and spiral galaxies with similar scatter, contrary to 
the A/bh-^b correlation (McLure & Dunlop 2002). Moreover, 
the SMBHs recently claimed to be detected in the globular clus- 
ters M15 and Gl, also appear to satisfy the same A/bh -o relation 
(Gebhardt et al. 2002). 



As both the w c - a and A/bh - a correlations seem to hold over 
the entire Hubble range (see e.g. de Zeeuw 2003 for a critical note), 
we can combine them to derive a correlation between the circular 
velocity and SMBH mass. For all galaxies in our sample, we have 
estimated A/bh via the A/bh —a correlation given by equation 
The results can be found in table 1. Combining our best fitting v c - 
a relation J5J with the A/bh —a relation J5J, we obtain 

(£)-(4* ± aWi 



lo. 



+ (7.24 ±0.17). (6) 



In Figure |2| (left panel) we plot the circular velocity versus the 
SMBH mass for the 24 galaxies of our sample with a > 80 km s ~ 1 : 
we find indeed a very tight correlation between A/bh and v c . More- 
over, this correlation also holds for the elliptical galaxies of Kron- 
awitter et al. (2000), as shown in the right panel of Figure|2| For the 
least massive galaxies with a < 80 kms -1 , this correlation again 
breaks down. 

The correlation between SMBH mass and circular velocity is 
useful for two different goals. Firstly, if a prescription can be found 
to link the circular velocity of a galaxy to the total dark halo mass 
A/dm, we can transform this Mbh-Wc correlation to a direct corre- 
lation between SMBH mass and total galaxy mass, which can then 
be compared with theoretical models of galaxy formation. A con- 
version between these two quantities can in principle be derived 
from high-resolution CDM cosmological simulations (e.g. Navarro 
& Steinmetz 2000; Bullock et al. 2001). Using the simulation re- 
sults of Bullock et al. (2001), 

3.32 

1.40 I — ) , (7) 



Mi 



DM 



10 12 M 



Ferrarese (2002) convert v c to A/dm, and derives a relation between 
SMBH and dark halo masses. If we repeat the same exercise with 
our larger sample of spirals and the most up-to-date A/bh - a cor- 
relation, we find 



A/ B 



0.11 



A/d 



(8) 



10 8 M \ ±ir- Mg ( 

It should be noted however, as indicated by Ferrarese (2002), that 
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the uncertainties in this conversion can be quite large. For exam- 
ple, a major uncertainty is how the baryonic infall in dark mat- 
ter haloes affects the baryonic rotation curve (e.g. Seljak 2002). 
Therefore, the correlation {5} shouls be regarded as a rough guide- 
line only. The A/bh - v c correlation JSJ on the other hand is based 
solely on observed quantities, and has a much smaller uncertainty. 
This correlation therefore serves as an important (purely observa- 
tional) constraint that should be reproduced and explained by the- 
oretical galaxy formation models. Combined with other tight rela- 
tions such as the A/bh - a relation and the Tully-Fisher relation, it 
clearly points at an intimate interplay between the various compo- 
nents (dark matter, discs, bulges and SMBHs), and forms a strong 
test for galaxy formation and evolution models. 

Apart from being an ingredient in theoretical galaxy formation 
models, the derived A/bh - v c relation can also serve as a practical 
tool to estimate the black hole content for spiral galaxies. Due to 
the large scatter of the A/bh - Lb relation, the most reliable means 
of estimating A/bh in galaxies is by using the A/bh - a relation. Un- 
fortunately, the number of spiral galaxies with reliable velocity dis- 
persion measurements is relatively small. On the contrary, extended 
rotation curves have been measured for large samples of spiral 
galaxies, mainly for use in Tully-Fisher relation studies (e.g. Persic 
& Salucci 1995; Courteau 1997; Verheijen & Sancisi 2001). This 
makes the tight A/bh - v c correlation a practical tool to estimate the 
black hole content of these galaxies. In particular, the A/bh -fc cor- 
relation can be used in a statistical way for SMBH demographics 
studies (e.g. Merritt & Ferrarese 2001c; Ferrarese 2002b; Aller & 
Richstone 2002; Yu & Tremaine 2002), by combining it with spiral 
galaxy velocity functions (Shimasaku 1993; Gonzalez et al. 2000). 
This falls beyond the scope of this Letter, and will be addressed in 
future work. 
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